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Abstract 

Recently, it has been pointed out that there exists a connection between the 
generalized parton distribution E(x, 0, t) and the Sivers asymmetry. For trans- 
versely polarized nucleon targets, generalized parton distributions are asym- 
metric in impact parameter space. This impact parameter space asymmetry, 
together with the final state interaction of the active quark, gives rise to the 
Sivers asymmetry in momentum space. We demonstrate this phenomenon ex- 
plicitly in the scalar diquark model. This result also illustrates the physics 
that underlies the correlation between the anomalous magnetic moment and 
the Sivers asymmetry for a given quark flavor. 
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1 Introduction 



Recently, one gluon exchange in the final state interactions (FSI) has been suggested 
[1] as a mechanism for generating a transverse single-spin asymmetry (SSA) in semi- 
inclusive deep inelastic scattering processes. This FSI can be effectively taken into 
account by introducing an appropriate Wilson line phase factor in the definition of the 
distribution functions of quarks in the nucleon [2, 3, 4, 5]. The physical interpretation 
of this Wilson line phase factor is that it describes the phase factor of the propagator 
for the active quark as it leaves the target. Therefore, the Wilson line extends from 
the position of the quark along a light-like future-oriented direction. This phase factor 
is not invariant under time- reversal, which is the reason why the Sivers asymmetry 
can be nonzero when the phase factor is included. 

This mechanism has a nice physical interpretation in transverse position space, 
where one can show that the Sivers asymmetry arises from a left-right (relative to 
the nucleon spin) asymmetry of the quark distribution in impact parameter space. 
The Wilson line phase factor describes the effect of the transverse component of the 
force that is acting on the active quark after it is knocked out of the target. This 
force is on average directed towards the center of the nucleon. This chromodynamic 
lensing effect thus translates the T-even transverse position space asymmetry into a 
T-odd transverse momentum space asymmetry of the leading quark [6, 7], which is 
the Sivers asymmetry. 

In Ref. [1] a simple scalar diquark model was used to demonstrate explicitly that 
the FSI can indeed give rise to a leading-twist transverse SSA, which emerged from 
interference between spin dependent amplitudes with different nucleon spin states. 
In Refs. [1, 8] it was observed that the same overlap integrals between light-cone 
wavefunctions that describe the anomalous magnetic moment contribution from a 
given quark flavor also appear in the Sivers distribution for that quark flavor (with 
additional pieces in the integrand). Since these integrals are the overlaps between 
light-cone wavefunctions whose orbital angular momenta differ by AL Z = ±1, the 
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orbital angular momentum of the quark inside the proton is essential for the exis- 
tence of the Sivers asymmetry. In Refs. [6, 7], the transverse distortion of impact 
parameter dependent parton distributions for transversely polarized targets in the 
infinite momentum frame, was used to develop a physical explanation for the sign 
of the Sivers asymmetry. In this paper, we diagonalize the relevant amplitudes in 
Ref. [1] by transforming to an impact parameter space basis [7]*. This should not 
only yield a clear intuitive physical interpretation for the SSA derived in Ref. [1], but 
at the same time provide an explicit example that illustrates a proposed connection 
between generalized parton distributions and SSAs [7]. 



2 GPDs in impact parameter space (scalar diquark 
model) 

Light-cone wavefunctions provide a very convenient representation of generalized par- 
ton distributions (GPDs) in terms of overlap integrals [13, 14]. For example, in a 
scalar diquark model in which a nucleon is modeled as a bound state of a quark and 
a scalar diquark, one finds in the case £ = (purely transverse momentum transfer 
A = p — p 1 ) 



H(x,0,-A 2 ± ) 

J_ f d2k 



(A x - iA v , 



2M 

1 r d 2 k 

4W (2?r 



(27T) 

E(x,0,-A 2 ± ) 



J (Z7T) 2 ^2 2 2 



, (1) 



, (2) 



where 

k^= k ± - (1 - x)A ± . (3) 

*For a general discussion of parton distributions in impact parameter space and their connection 
to generalized parton distributions the reader is refered to Refs. [9, 10, 11, 12]. 
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For higher Fock components, similar convolution formulas exist [13, 14]. 

A well known feature of Fourier transforms is that they diagonalize convolution 
integrals. For this purpose we switch to the transverse position space representation 
of the light-cone wavefunction 

d 2 k 



i) s (x, Cj.) = 



(2^ 



Jk ± c ± 7\ 



(4) 



Note that the transverse momentum in the two-particle Fock component is Fourier 
conjugate to the distance Cj_ = rj_i — r± 2 between the active quark and the spectator. 
However, since GPDs have a particularly simple form in impact parameter representa- 
tion, we change variables from to the impact parameter (the distance between 
the active quark and the center of longitudinal momentum), which are related by 



(1 - X)C L , 



(5) 



yielding 
H(x,b ± ) 



d 2 A 



(27T) 



1 

47T 



Likewise one finds 

„• d a 



' db x dbv 

2M 



£{x,b±) 



1 

47T 



^lt(x, cj.)^i.i(a:,cj.) + iplt(x, cj.)^ii (x, c ± ) 

"To "To o o 



(6) 



1 — x 



where 



S(x,b ± ) = J 



d 2 A 



(2tt; 



L e iA ± -bj 



E(x,0,-A 2 ± ) . 



(8) 



The physical significance of £(x,b±), i.e., the spin-flip distribution in impact 
parameter space becomes clear when we consider a state that is polarized in the +y 
direction (in the infinite momentum frame) 



P + , R± = 0±, +y) = ^= (\P + , R± = ± , t) + % \P + , R± = 0±, |)) . 



(9) 



For this state, the unpolarized quark distribution in impact parameter space reads 



q y (x,b±) 



P + ,R± = 0±,+y 
d 2 A 



(27T) 



JA. ± -b ± 



O q (x,b ± ) P+,R ± = ± ,+y 
H(x, 0, -A 2 ± ) + i^— E(x, 0, -A 



2M 



H(x, b±) + 



1 d 



£(x,b ± ) . 



(10) 



2M db x 

We will write in Section 4 on the interesting (and important) property that the 
density in impact parameter space qy(x, h±) has the asymmetry along the ^-direction 
in impact parameter space when the proton spin is polarized along the y-direction 
as we see in (10), whereas there does not exist such an asymmetry for the density in 
momentum space. 



3 Transverse distortion of the wavefunction 

In Ref. [12] it was shown that if the spin-flip generalized parton distribution E q (x, 0, — A 2 ± ) 
is nonzero, then the parton distribution of quarks with flavor q is distorted in the 
transverse plane when the target has a transverse polarization. For a nucleon with 
spin pointing in the positive y-direction, one finds [12] 

1 d 

q y (x,b ± )=H(x,h ± ) + ——£(x,b ± ), (11) 
where Ti and £ are Fourier transforms of generalized parton distributions 
n q (x,b ± ) ee J^e^ b -H q (x,0,-Al) 

S 9 (x,b ± ) = J^e^ b -E q (x,0,-Al). (12) 

One of the things that are known about £ q is that its integral should give the contri- 
bution from flavor q to the anomalous magnetic moment^ 

J dx J d 2 b L £ q {x,b ± ) = Kq . (13) 

^The measured magnetic moment of the proton is obtained as a superposition from all quark 
flavors Hp — J2 q e q K q = 1-79. 



Integrating E(x, 0, — Aj_) over x yields the Pauli form factor F 2 (— Aj_). Therefore, 
unless the x dependence introduces a fluctuating sign into E(x, 0, — Aj_), the sign of 
F2 thus determines the sign of E(x, 0, — Aj_). Hence one expects that E(x, 0, — A^_) 
is a smooth function of Ai with a maximum (or minimum if K q < 0) at the origin. 
From experience with Fourier transforms it is thus clear that £(x, h±) also has a 
maximum (minimum) at the origin and is otherwise a smooth function that will have 
the same sign as n q for all (or most) values of bj_. 

The fr^-derivative of a smooth positive function S(x, h±) with maximum at the 
origin (we consider here the case K q > 0) is positive for negative b x and negative for 
positive b x . For valence quarks, 7i(x, h±) is known to be positive. Therefore, adding 
the b x - derivative of £(x,h±) to H(x,h±) has the effect of shifting the distribution 
towards negative b x . For negative K q the effect is reversed. This is why the sign of n q 
determines the sign of the distortion of the quark distribution in impact parameter 
space. For a nucleon that is polarized in the +y direction the distortion is towards 
negative x when K q > and towards positive x when n q < 0. 

Note that although the argument hinges somewhat on assumptions about the 
shape of E q {x, 0, Aj_), these assumptions seem to be satisfied for typical model ansatze 
for GPDs and therefore the result is actually rather general. In the specific example 
of the scalar diquark model, we can calculate £ q (x, b±) and verify that the model 
satisfied the above assumptions. 

Since the scalar diquark model provides us also with the light-cone wavefunction 
for the 'nucleon' there is actually a more direct way to determine the sign of the 
distortion in impact parameter space. For this purpose we consider the wavefunction 
of a quark that is polarized in the +y-direction [15, 16] 



^i(x,k ± ) 



2 




(14) 



(15) 
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1 

71 



i[M + 



m 



x 



k x + ik y 



x 



4>(x,k ± ). 



Although this should be evident from time-reversal invariance, we note that the naive 
(i.e. gauge noninvariant) unintegrated momentum space distribution obtained from 
the wavefunction [described by Eqs. (14,15)] squared is even in k^. Indeed 



qy(x, kj 



4ttV 



V3 ( x ' kj 

2 



+ 



V3(x,kj 

2 



1 

47T 



mV k, 
M+- +^ 

x / x^ 



2 , (16) 



where we used the fact that is real. (See Eq. (38) in Appendix A.) 

However, a transverse asymmetry in position space is not excluded by time- 
reversal invariance and Eqs. (14,15) do in fact correspond to a state with an asym- 
metry in the x-direction as we will now demonstrate explicitly. After performing a 
Fourier transformation to the transverse relative position space coordinate cj_, we 
have 



with 



tp+\(x,c±) = J 



2 



i 

d 2 k 



ic i -k 



i/3 ( x , k J 

2 



m Id id 

x x dc x x dc y 



0(C±) 



(2tt 



e ic± ' kl 



^3 ( x > k J 



1 

71 



id Id , ( , , m\ 

x ac^ x ac^ V x / 



0(c±) 



(17) 



(18) 



0(c ± ) = 



d 2 k ± 
(2^ 



3 ic ± -kj_ 



= -gxy/l - x J 



0(k±) 
rf 2 kj 



ic ± -k i 



k 2 ; +5 



-^f ^v / r^^o(|c±iv / 5) , 

Z7T 



(19) 



and 



B 



9 m 2 A 2 . 
= x 1 - x -M + — + 

x 1 — x 



(20) 



In the limit \c±\y/B — > 0, (19) becomes (g/^Tr) xy/1 — x ln(c\B). Using (17) and 
(18), we have 



4>H(x, Cj.) 

2 



with 



ipVi(x, c ± ) 

2 



^+|(x,c_l; 

2 



1 

2 



^3(x,c ± ) 



1 d 



) 



(21) 



M+-U+-— , 

x J x dc x 



+ 1 ( d 

x 2 1 dc 1 



+ - M+-U— 0, 



(22) 

since 0(c_|_) is real. In (22) the last term is odd under c x — > — c x and describes the 
deformation of the target in impact parameter space as predicted by Eq. (11): 



1 9 e( u ^ 1 



2Mdb x 



47T (1 — X) 2 X 



M 



m 



X 



> d 

dc 



:0- 



(23) 



Since <j>(x, c±) is a monotonically decreasing function of c 2 ± , this implies that 



ip+l(x,c ± ] 



distorted towards negative c x . For example, the mean _L coordinate for 



fixed x yields 



<c*> = / 



d 2 c L 

1 1 



47T 



47r 2x V rr 



c 

'2„ jl2 



(24) 



cf c ± z < . 



4 SSA in impact parameter space (scalar diquark 
model) 

The SSA for the scalar diquark model has been calculated in Ref. [1]. For a target 
that is polarized in the +y direction, one finds for the _L momentum distribution of 
the outgoing quarks* of quarks that carried light-cone momentum fraction x 

VH^±) /i(x,rj.) 

*Note that there should be an additional overall minus sign in front of Eq. (21) of Ref. [1]. 
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= c "^"'"-"^/|^ (k , ( - "ff v - (25) 

where C = g 2 eie 2 /(2(27r) 3 ) and j3 is given in (20). Here P^(x, rj_) is the actual spin 
asymmetry and fi(x,r±) is the unpolarized quark distribution so that the product 
V^fi is the spin-odd part of the outgoing quark momentum distribution. 

In the following we are going to investigate within the context of the scalar diquark 
model if this SSA can be related to GPDs and the asymmetry of PDFs in impact 
parameter space, as has been conjectured in Ref. [7]. Eq. (25) yields for the average 
transverse momentum in the x direction 

(r x ) = j d 2 r± r x V § (x,r ± ) fi(x,r ± ) (26) 

= c (I « + m) (1 - ,) / / *l ^ ^ p^fe 

Using Eq. (19) 

kiTB-WT^/^'" 1 ^* (27) 



2 

9 



and 

1 



±== J d 2 c ± e-^ c ^(c±) (28) 



r 2 + B gxy/1 - x 
as well as the relation (46) given in Appendix B, from (25) we get 

ei e 2 (xM + m) r 2 / Ai f ,o „_ ir± . cj /-,2/ ik,-c 



V 7 2 2vr 3 x 2 7 7 2tt 2 7 7 



2(2tt) 3 x 2 7 7 (2tt) 

9\, , (fc±-r ± ) 



v 9cx/ (k ± -r ± ) 2 + A 9 

-ieie 2 (xM + m) f 2 . 9 . r d 2 k ± - k fc^ 



x 0(c'_l) ( -i— ) 0(c± y/1 . . x 2 

2(2,) * / A ^»^»/f^ e ' kiCi ki + A/ 

/ rf 2 b ± I±(c±)— ^-5(x,b ± ), (29) 
7 ±v 2M db x v ; v ; 



where 



^> - -^/f^- ^ (30, 



a 9 ^o eie 2 d 



4tt c 2 . 



(31) 



is the transverse impulse as a function of the _L coordinate of the struck quark. 

The physical interpretation as a 'transverse impulse' becomes clear when one 

compares Ij_ to the net transverse impulse / dtF± that one obtains when one integrates 

the transverse component of the Coulomb force along a straight line from position 

(c^,Cj_,0) to (c^,Cj_,+oo) along a straight line (remember: for an ultrarelativistic 

particle, z(t) = t) 

f°° 7, m / ,,h eie 2 [°° c ± e ie2 c ± 

/ dt F ± (c ± ,z(t)) = — — / dt r = — T (32) 

Jo 4tt Jo ( c ^ + t 2)l 4tt c\ v ' 

We note that Eq. (29) can be written as 

(O = /^Vil(c)^A^ b±) 

= J d% ± F|(c ± ) [qf™(x, b ± ) + q^(x, b ± ) 

= Jd 2 b ± V ± (c ± )q y (x,b ± ), (33) 

where q y ym {x, bj_) and q^ yra {x, bj_) are symmetric and asymmetric about b\ = 0, and 
they are given by qf m (x,b ± ) = H(x,b ± ) and qf yni (x, b±) = ^-^£(x,b±). This 
relation between GPDs and SSAs, which we now derived explicitly for the scalar 
diquark model, is exactly of the form U SSA = GPD * FSF that was proposed in 
Ref. [7]. 

We emphasize that q^ 7m (x, bj_) is not zero in the impact parameter space and this 
fact made the expressions given in Eq. (29) and (33) possible, however, q^ sym {x, kj_) 
in the momentum space is identically zero: 

qy(x^±) = qj m (x,k ± ) , (34) 

as we demonstrated explicitly in Eq. (16) in the scalar diquark model which has the 
wavefunctions given in Eqs. (37) and (39). Even though there exists such a difference 
between the asymmetric parts of q y (x, bj_) and q y (x, kj_), they satisfy the following 
relation which is consequence of a general property of Fourier transforms: 

d 2 b ± q y (x, b±) = J q y (x, kj_) , (35) 
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i.e., the norm (integral of the absolute square of the function) is invariant under a 
Fourier transform. 

Our calculations also illustrate explicitely a very important difference between 
impact parameter dependent parton distributions q y (x, b±) and unintegrated parton 
distributions q y (x, kj_) for a transversely polarized target: Even for a transversely 
polarized target the unintegrated parton distribution q y (x, kj_) has to be symmetric 
under kj_ — > — kj_ because of parity and time-reversal arguments. Roughly speaking, 
the reason is that given only the spin vector S and the momentum vector P of the 
target it is not possible to construct a parity and time-reversal invariant correlation 

— * — » — » 

between the target spin and the quark momentum k that is odd under k — > —k. For 
example, although the term k - (^S x P^j is invariant under parity, it is odd under time- 
reversal and therefore cannot appear in the momentum distribution for the quark§. 
However, the situation is different for the transverse position bj_ where an asymmetry 
with respect to — > — bj_ is allowed: To see this we consider the term b ■ (^S x P*j 
which is clearly invariant under both parity and time-reversal. The fact that one can 
write down a term proportional to the spin that is odd in b and consistent with the 
symmetries of QCD means that a transverse (relative to P and S) asymmetry in the 
position space distribution is allowed. 

The explicit model calculations in the scalar diquark model confirm these general 
considerations since they yield q y (x, — kj_) = q y (x, kj_) [Eq. (16)], while q y (x, — bj_) ^ 
q y (x,b ± ) [Eq. (24)]. 

We should also point out that Eq. (33) is consistent with expressions that have 
been written down previously for the mean transverse momentum arising from the fi- 
nal state interactions [17, 18]. What is new is that we are evaluating these expressions 
in impact parameter space and for a specific model for the final state interactions. 



§ This argument is also the reason why a (time-reversal invariance breaking) initial or final state 
interaction is crucial for the Sivers asymmetry in QCD. 
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5 Summary 



We have analyzed the single-spin asymmetry in semi-inclusive deep inelastic scattering 
(DIS) in the scalar diquark model. In the impact parameter representation, the SSA 
emerges as a correlation between the distribution of partons in the transverse plane 
and the transverse impulse, which the quark being ejected from a certain transverse 
position has acquired as a result of the final state interactions. 

The scalar diquark model is a model where both the light-cone wavefunctions 
of the quarks as well as the final state interactions are constructed perturbatively 
and therefore the constructed amplitudes still have all Lorentz symmetries as well as 
gauge invariance built in. Because of all these features, the physics of the Wilson line 
phase factor that describes the final state interaction in semi-inclusive DIS should 
be correctly represented in this model. In this model we obtained in the impact 
parameter space the formulas (29) and (33) which show explicitly the connection 
between the SSA and £(x, h±) which is the Fourier transform of the generalized parton 
distribution E q (x, 0, — Aj_). These formulas are also useful for understanding the role 
of the orbital angular momentum in the Sivers asymmetry, since E q (x, 0, — A^_) is 
given by the overlap integrals between light-cone wavefunctions whose orbital angular 
momenta differ by AL Z = ±1 [1, 14, 16]. 

Without the final state interaction, one should not find a transverse momentum 
asymmetry because of time-reversal invariance. However, time-reversal invariance 
does not exclude a transverse position space asymmetry when the nucleon is polar- 
ized transversely. Our explicit calculations within the context of the scalar diquark 
model confirm these general predictions. We also find that the transverse position 
space asymmetry is described by the Fourier transform of the generalized parton 
distribution E(x, 0, — A^_) as one would expect [12]. The transverse single-spin asym- 
metry in the model is obtained by convoluting the final state interaction kernel with 
the momentum space light-cone wavefunctions of the 'quarks' in the 'nucleon'. Upon 
Fourier transforming the SSA to impact parameter space, one finds that the average 
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SSA has a probabilistic interpretation in impact parameter space in the sense that 
the asymmetry can be obtained by averaging the transverse impulse for each point in 
impact parameter space with the probability density to find a quark at that impact 
parameter. On the one hand, this result illustrates clearly the physical mechanism 
for the results in Ref. [1], and on the other hand it provides a specific example that 
confirms the general mechanism for SSA, that was proposed in Refs. [6, 7]. 



One of the virtues of light-cone wavefunctions is that many form factors and transition 
matrix elements have very simple representations as overlap integrals involving these 
wavefunctions [19, 20]. Of course, in order to apply these overlap integrals, one needs 
to know the light-cone wavefunctions for each Fock component. Since we do not 
know the light-cone wavefunctions for the nucleon in QCD, but would still like to 
use them to illustrate the connection between transverse single-spin asymmetries and 
impact parameter space distributions, we take them from a simple toy model (a scalar 
diquark model). 

In the scalar diquark model, the light-cone wavefunction for the two-particle Fock 
component is obtained by calculating the 'splitting' of a fermion into another fermion 
plus a scalar perturbatively. This yields the two-particle Fock component wavefunc- 
tions of the J z = +\ state [15, 16]: 



A Scalar diquark model 




where 




(37) 
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and 

* = & x > k ^ = 7!= M , g±g El±g ' (38) 

In general one normalizes the Fock state to unit probability. 

Similarly, the J z = — \ two-particle Fock state has components: 



^ii(x,k ± ) = (M+^)0 . 



(39) 



The spin-flip amplitudes in (37) and (39) have orbital angular momentum projec- 
tion l z = +1 and —1 respectively. The dependence in the numerators of the 
wavefunctions is characteristic of the orbital angular momentum, and holds for both 
perturbative and non-perturbative couplings. 



B Useful relations for Fourier transformations 

The fact that certain amplitudes that are convolutions in momentum space become 
diagonal in position space can be easily understood on the basis of some elementary 
theorems about convolutions and Fourier transforms. For example, if 

/(k ± ) = J<Pb ± e-*^f(b ± ) 

g(k ± ) = | A i e- ik - b ^(b 1 ) (40) 

then the "form factor" 

F{q±) = I |^ r(k±) ^ (k± + q±) (41) 
becomes diagonal in Fourier space 

/ |^ 6iq± ' b±F(q±) = ^ b ^(bJ. (42) 

This well-known result forms the basis for the interpretation of non-relativistic form 
factors as charge distributions in position space. 
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Similarly, the convolution describing the FSI in the amplitude has the form 

A(k ± ) = J !^(^ - q^(qx) (43) 
and thus becomes diagonal in position space 

where 

K(b ± ) = J |^(<uK* q± - bx (45) 

is the Fourier transform of the kernel. 

Finally, the mean SSA in momentum space contains a double convolution: 

= /^bx/'tbxj^fbxjsfbx). (46) 

Of course, for practical purposes it is often more convenient to stay in momentum 
space since experiments measure momenta and Fourier transforms are numerically 
awkward. However, for the purpose of physical interpretation, the position space 
expressions are useful since they are diagonal and therefore allow a probabilistic in- 
terpretation. 
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